Quantum spin Hall effect induced by non-magnetic and magnetic staggered potentials 
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We have a comparative study of the quantum spin Hall (QSH) effects induced by non-magnetic 
and magnetic staggered potentials respectively and show that they have the same effect in driving 
the topological phase transition. The result implies that both time-reversal (T) preserving and 
breaking systems can host QSH effect. We also investigate the stability of the resulting QSH effect 
to disorder and find that for T invariant system the edge states are always robust while those of T 
breaking system are also robust if there is additional symmetry in the system. 



PACS numbers: 73.43.-f, 72.25.Hg, 73.20.-r, 85.75.-d 

Recently the field of topological insulator (TI) has at- 
tracted a great deal of interests, due to their exotic phys- 
ical properties as well as potential applications, such 
as spintronics, quantum computing, et al PHI]- Many 
materials have been predicted and discovered to show 
TI phases (including HgTe/CdTe quantum wells (QWs) 
[HE], bismuth antimony alloys [3 [8]; Bi 2 Ses, Bi 2 Te3, 
and Sb 2 Te 3 [5HT2]; Heusler compounds [H[H]; Tl-based 
ternary chalcogenide series [T5J [T5] ; et al) . The findings 
of the real materials not only provide a platform to test 
the predictions of many unusual phenomena exhibited by 
TI [17H20], but also inspire more theoretical studies on 
TI. 

The study of TI begins at the theoretical proposal of 
QSH effect in graphene by Kane and Mele [HI 122 ■ They 
expected that spin-orbit coupling will convert graphene 
from an ideal 2D semimetallic state to a QSH insula- 
tor. The resulting QSH insulator is topologically distinct 
from a band insulator, so it is referred as TI. However the 
calculations have suggested that the spin-orbit coupling 
in graphene is too small to reveal the QSH effect exper- 
imentally [23] • Remarkably, in 2007 the QSH effect was 
realized in HgTe/CdTe QWs following the theoretical 
suggestion of Bernevig, Hughes and Zhang [5] [6] . Later 
though many studies have been carried out in identify- 
ing new physical systems that will possess topological 
nontrivial phases, till now the 2D TI is only found exper- 
imentally in HgTe/CdTe QWs. In the low-energy effec- 
tive theory, the QSH effect can be understood from Dirac 
Hamiltonian with masses [24j [25] . The relative signs of 
the masses at the Dirac points determine the phases of 
the system. Alternately, it can also be understood from 
band inversion, which is the mechanism of the TI phase 
in HgTe/CdTe QWs. The two ways are equivalent since 
the occurrence of band inversion corresponds to changing 
the sign of one Dirac mass and causes a topological phase 
transition, which can't happen without closing the gap. 

Spin-orbit coupling is a necessary condition for the ex- 
istence of TI. Its role is to induce a gap in the Dirac dis- 



persion and to ensure that the gap is finite everywhere in 
the Brillouin zone (BZ). It also has been known that non- 
magnetic and magnetic staggered potentials can perturb 
the Dirac dispersion and induce a gap. Such terms can 
be obtained by the proximity effect to the corresponding 
orders (charge density wave (CDW) and antiferromag- 
netism (AF)) [26]. When these terms coexist, their in- 
terplay will determine the phase of the system. In this 
paper, we study the interplay of these terms. We start 
from a trivial insulator with spin-orbit coupling and in- 
troduce non-magnetic and magnetic staggered potentials 
with check-board and stripe patterns into the system. 
We find that when the strength of the potential is strong 
enough a band inversion occurs and the system shows 
QSH effect. Specially the QSH effect induced by mag- 
netic staggered potential breaks T symmetry, in contrast 
to the previously studied QSH effect. 

To be concrete, we study a model describing 
HgTe/CdTe QWs. It resides on a square lattice with 
four orbit states |s, f), \Px + iPyA), \s,i), \ {p x ~ iPy),i) 
(f, I denote the electron's spin) on each site. In the mo- 
mentum space, the Hamiltonian writes, 

iZo(k) = [4D-2D(cosfc :E +cosfc y )]/ (1) 
+ [M + AB - 2B(cos k x + cos k y )]a z 
+ 2 A sin k x s z (8 <j x + 2A sin k y cr v 

Here? and s are Pauli matrices representing the or- 
bits and the electron's spin and I is identity matrix. 
A, B, D and M are four independent parameters. The 
tight-binding Hamiltonian can be directly obtained by 
a lattice regulation of the effective low-energy Hamilto- 
nian describing the physics of HgTe/CdTe QWs. We 
can also view it as a simple toy model conveniently de- 
scribing both topological and ordinary phases of non- 
interacting electrons in 2D. The energy spectrum of 
Ho(k) has two double degenerate branches = (AD — 

D k ) ± y '(2Asmk x ) 2 + (2A sin k v ) 2 + (M - B k ) 2 , where 

M = (M + AB), S k = 2B(cosk x + cosk y ) and D k = 
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FIG. 1: (Color online). The staggered potential with check- 
board pattern (a) and stripe pattern (c) and the band evolve- 
ments at TRIM ((b) and (d)) with the potential strength cor- 
responding to the patterns in (a) and (c) respectively. The 
open and filled circles on lattice sites in (a) and (c) represent 
on-site potentials with equal values but opposite signs. The 
new lattice vectors and new BZs are shown in (a), (c) and the 
insets of (b), (d). The blue curve in (d) is a band evolving at 
other momentum which will determine the gap size at large 
potential strength. The range with vertical dashed lines in 
(b) and (d) marks the gap of the induced topological phase. 
The parameters are fixed for all calculations in this paper 
to be A = 36A5meV, B = 27.44meV, D = 20.48meV and 
M = ImeV. 



2D(cosk x + cosky). At half-filling, depending on the 
values of M and B, the system can be QSH or trivial 
insulator. 

The system defined by Eq. ([I]) is invariant under 
T and spatial inversion. Since the inversion operator 
V operating on p— type orbit generates a minus sign, 
the inversion operator writes V — I ® a z . Thus the 
signs of M — _Bk at the four time-reversal invariant mo- 
menta (TRIM) determine the phase of the system, i.e., 
M(M + 4B) 2 (M + 8B) < 0(> 0) for QSH (trivial) insu- 
lator [221 [21] ■ In this paper, we restrict our calculations 
in the parameter range describing HgTe/CdTe QWs, 
where B > \M\ > 0. Since the TRIM T 1 (k = [0,0] 
in the BZ, see the inset of (b) in Fig. ([I])) dominates 
the physics in the range, the Hamiltonian describes QSH 
effect for M < and a trivial insulator for M > 0. Ex- 
perimentally the gap parameter M can be continuously 
tuned from a positive value for thin QWs with thickness 
d < d c to a negative value for thick QWs with d > dd (d c 
is a critical thickness and equals 6.3nm for HgTe/CdTe 
QWs) [B]. The gap parameter M in Eq. |l]) repre- 



sents an on-site potential, which has different sign for 
s— type and p— type orbits. When changing its sign, one 
of the occupied bands changes from p— type (M > 0) to 
s— type (M < 0). Since the two kind of orbits have dif- 
ferent parities, a band inversion will induce a topological 
phase transition. 

To drive the system into QSH phase, a band inversion 
is needed. So it is interesting to seek ways other than 
tuning the gap parameter M to generate the band inver- 
sion. Below we fix M > when the system is a trivial 
insulator and find ways to inverse the bands at TRIM 
T\. A natural thought is to enlarge the Hamiltonian, 
which make it possible to add more terms to it. Firstly 
we consider putting the system on a check-board square 
lattice (CSL), which can introduce alternating potential 
with check-board pattern. In the case, the unit-cell is 
doubled. The Hamiltonian -ffo(k) is enlarged to 8 x 8 
and becomes, 

Hi (k) = ADI - 2£>(cos k x + cos k y )t x (2) 
+(M + 4B)<r z — 2H(cos k x + cos k y )r x ® a z 
+2A sin k x r x ® s z <E) o x + 2 A sin k y r x ® a y 

Here f is Pauli matrix describing the two sublattices. We 
have been able to identify two interesting on-site terms: 

(1) a non-magnetic staggered potential (or CDW poten- 
tial) V\T Z \ (ii) a magnetic staggered potential (or AF po- 
tential) V%t z <8) s z (The magnetization can also lie in the 
plane and couples with the in-plane components of the 
electron's spin, generating terms like t z <g) s x or t z ® s y . 
We will discuss these terms later) . The former term pre- 
serves T symmetry while the latter breaks. Including the 
above terms to Eq. though we can't obtain the an- 
alytic forms of the energy spectrums, they are the same 
for the two different cases. The system remains gapped 
and the occupied bands evolve with the strength of Vj 
(i = l,2). 

At TRIM r l5 the energy eigenvalues are: (1) 4D-M- 
(2) 4D- M+B-; (3)AD+M -B+; (4) 4D + M+B+, 
where B± = a/(4D ± 4B) 2 + V? and each is double de- 
generate. The corresponding eigenvectors and their pari- 
ties can also be obtained, which are listed in Table I. The 
band evolvement with the strength Vi is shown in (b) of 
Fig. 0. At half-filling and for Vi = 0, bands (1) and 

(2) are occupied. Since we choose M > 0, the system 
is a trivial insulator. As Vi increases, bands (2) and (3) 
firstly approach each other and at a critical value of Vi 
the filling for the two bands will interchange. Bands (2) 
and (3) consist of electrons in the p— type and s— type or- 
bits respectively and have opposite parities. For the case 
with CDW term the band inversion will induce a topo- 
logical phase transition and drive the system into QSH 
phase. It is interesting that AF term has the same effect 
as CDW term. As we further increase the strength of the 
potential, the system remains in QSH phase until another 
band inversion occurs at other TRIM, which happens at 
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No. 


Eigcnvector(CDW) 


Eigenvector(AF) 


Parity 


1 


(-0l_,l 6 ) 


(-^-,i 8 ) 


(-0L,i 6 ) 


(-0t-,i 8 ) 


-1 


2 


(0?-,l 6 ) 


(^-,i 8 ) 


(0?-,i 6 ) 


(^-,i 8 ) 


-1 


3 


(0i+,i 7 ) 


(02+! I 5 ) 


(0?+,i 7 ) 


(02+, l 5 ) 


1 


4 


(-0?+,i 7 ) 


(-0i+,i 5 ) 


(-0i+,i 7 ) 


(-01+, I 5 ) 


1 



TABLE I: The eigenvectors for Hamiltonian at TRIM I\ and 

their parities. Here cj>i± — ±^±41) an( ^ 2 ± ~ 4B±4D ■ The 
superscripts on the values of the eigenvectors represent the 
corresponding position in the eigenvectors. 



potential strength bigger than 100 meV. The gap range 
of the resulting QSH insulator is also denoted in (b) of 
Fig. ([lj. At small potential strength, it is determined by 
band (2) and (3) and increases with increasing potential 
strength. Then a band evolving at TRIM T 3 (k = [0, tt] 
in the BZ, see the inset of (b) in Fig. 0) becomes the 
lower restriction of the gap and restricts its further in- 
crease. It is interesting to note that the behavior is sim- 
ilar to what happens in topological Anderson insulator, 
where the phase diagram has to be obtained by conduc- 
tivity calculations since the system has no translation 
symmetry in the presence of disorder [29H32] . 

We can also put the system on a stripe square lattice 
(SSL), as shown in (c) of Fig. 0. The Hamiltonian 
Ho(k) becomes, 

i? 2 (k) = {AD - 2D cos k y )I - 2D cos k x r x (3) 
+(M + AB-2B cos k v )a z - 2B cos k x r x ®cj z + 
2 A sin k x r x ® s z ® a x + 2 A sin k y a y 

Similarly, a stripe CDW (r z ) or AF term (t z ®s z ), which 
has the same form as its counterpart on CSL, can be 
added to the above Hamiltonian. The energy spectrums 
for both cases are still the same. The band evolvements 
with the potential strength at TRIM are shown in (d) 
of Fig. 0. The CDW or AF term on SSL can also 
induce a band inversion and drive the system into QSH 
phase. However compared to the case in CSL, the band 
inversion here occurs at smaller potential strength and 
the resulting QSH insulator has bigger gap. 

So adding CDW or AF term to the system can induce 
a band inversion. Though they have different properties 
under T transformation, they have the same effects on 
band inversion. The reason is that the combined Hamil- 
tonian at TRIM Ti can be decoupled for each orbit and 
spin. Each decoupled Hamiltonian has a 2 x 2 form and 
has two sub-bands. The CDW or AF term pushes one 
sub-band up and the other down, which doesn't depend 
on the sign of the potential strength. Since the differ- 
ence between CDW and AF terms is that AF term acts 
on different spin with different sign, they have the same 
effect on changing the band structure. 

The CDW term preserves T symmetry and such sys- 
tems can be described by Z^ topological invariant [2"2"ll33| . 



Since our system has inversion symmetry, the Z^ invari- 
ant can be determined from the knowledge of the pari- 
ties of the occupied band eigenstates at the four TRIM. 
When the band inversion occurs as we increase the po- 
tential strength, the value of the Zi invariant also change 
its sign and becomes non-trivial, indicating that the sys- 
tem is in QSH phase. For the case with AF term, T 
symmetry is broken and Zi topological invariant is in- 
applicable. However the spin Hall conductance (SHC) 
of the resulting QSH phase still shows quantized value 
^ [51] . So the underlying topological invariant is a spin 
Chern number, which describes the quantized spin-Hall 
conductivity. The concept of spin Chern number has ap- 
peared in the recent literature which has its definition for 
systems with spin s z conservation [35] - For T invariant 
systems, it is equivalent to Zi topological invariant [33] . 

To further support our identification of the topological 
phase, we have performed numerical diagonalization of 
Hamiltonian (2) and (3) with CDW or AF term using a 
strip geometry in the range of parameters where the sys- 
tem is in band-inverted phase. In accord with the above 
arguments we find a pair of spin-filtered gapless states 
associated with each edge traversing the gap, which is 
shown in Fig. The spin-filtered edge states deter- 

mine the transport of charge and spin in the gap range. 
For a two-terminal device, the conductance is contributed 
by two conducting channels on the edges and gets quan- 
tized value 2e 2 /h. For a four-terminal device with proper 
voltage on each terminal, a spin current can be generated 

Up to now, we demonstrate the existence of QSH ef- 
fect induced by non-magnetic and magnetic staggered 
potential in a trivial insulator with spin-orbit coupling. 
In the following, we study the stability of the resulting 
QSH effect to non-magnetic disorder. The Hamiltonian 
we are considering can be decoupled for spin-up and - 
down electrons and each describes quantum anomalous 
Hall (QAH) effect. The resulting QSH effect can be un- 
derstood as two copies of QAH effect. For the CDW case, 
the two copies for spin-up and -down electrons are related 
by T and the QSH effect is immune to non-magnetic dis- 
order. However for the AF case, though T symmetry 
is broken, the system preserves the combined symmetry 
of T and a primitive lattice translation, which has been 
studied in three dimensions and antiferromagnetic topo- 
logical insulator is predicted [3B]. So the two copies for 
spin-up and -down electrons in the presence of AF term 
are related by the combined transformation. If there is 
non-magnetic disorder in the system, the combined sym- 
metry will be broken and the two copies will behave sep- 
arately. But we still expect the combined system will be 
robust to disorder, because each copy is in QAH phase 
and robust to disorder. 

To support the above statement, we employ the re- 
cursive Green's function method to evaluate the conduc- 
tance G of two-terminal devices (see the insets of Fig. 
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FIG. 2: (Color online) One-dimensional energy bands with a 
strip geometry (shown in inset) for (a) check-board pattern 
and (b), (c), (d) stripe pattern. In (a) and (b), the energy 
bands are the same for cases with CDW and AF terms, (c) 
and (d) have different edges with (b) and the energy bands 
with CDW term (c) are different with those with AF term 
(d). A strip of width N y = 60 unit cell with open bound- 
ary conditions along y and infinite along x is used with the 
staggered potential strength 30meV, when the system is in 
band-inverted phase. 



([2])) using Landauer-Biittiker formalism. Figure [3] shows 
the results of such calculations in the space of parameters 
(Ep, Uq), where Uo is disorder strength and the disorder 
is described by a random on-site potential uniformly dis- 
tributed in the range (—Uq/2, Uq/2). We only consider a 
single disorder realization at each point of the (Ep, Uo) 
phase diagram. Nevertheless, this turns out to be suffi- 
cient for studying the stability of the edge states. The 
reason is that if the edge states are robust to disorder 
there should be a region showing quantized conductance 
G = 2e 2 /h in the phase diagram and conductance G in 
the region shows no observable fluctuations, but fluctu- 
ates significantly elsewhere. The plots in Fig. ^ dis- 
plays conductance G in a fashion that is designed to am- 
plify the effect of fluctuations. In (a), (b), (d) and (e) of 
Fig. (|3| , regions showing no observable fluctuations exist 
in the phase diagrams, implying that the QSH phase in- 
duced by CDW or AF (in s z channel) terms is robust to 
disorder. We also carried out calculations on strip SSL2 
and the results are similar to those on strip SSL1. These 
results are consistent with what we expect. 

As mentioned, for AF term, the magnetization can also 
lie in the plane (generating term in s x or s y channel) 
and in the following we will focus on such term. Includ- 
ing such term to Hamiltonian described by Eq. (2) or 
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FIG. 3: (Color online) Conductance G as a function of dis- 
order strength Uo and the Fermi level Ef- Each data point 
corresponds to a single disorder realization, (a), (b) and (c) 
are the results with CDW term, AF terms in s z and s x chan- 
nels on CSL strip; (d), (e) and (f) are corresponding results 
on SSL1 strip. The parameters here are the same as those in 
Fig. pj, when the system shows band-inverted phase in its 
clean form. 




-71/2 k jt/2-j[H k nil -it -Till k it/2 71 



FIG. 4: (Color online) One-dimensional energy bands with 
AF term in s x channel on a strip geometry (a) CSL; (b) SSL1; 
(c) SSL2. The parameters used here are the same as those in 
Fig. (§. 

Eq. (3), the band evolvement at TRIM Ti is exactly 
the same as that due to CDW or AF potential in s z 
channel. However unlike the cases with s z channel term, 
the combined Hamitonian can't be decoupled for spin- 
up and spin-down electrons (i.e. there is no longer s z 
conservation in the system). We perform numerical di- 
agonalization of Hamiltonian (2) and (3) with AF term in 
s x channel (it is the same with s y channel term) and the 
one-dimensional energy bands with a strip geometry are 
shown in Fig. Q. For strips CSL and SSL1 (see insets 
of Fig. ([2]) ) , the energy spectrums are similar to those in 
Fig. ([2]) and there are edge states traversing the gap in 
band-inverted phase. But for strip SSL2, the edge states 
vanish though the system is in band-inverted phase. The 
difference can be understood from the symmetries exist- 
ing in the system. AF term in s x channel couples spin-up 
and -down electrons and the system only preserves the 
combined symmetry of T and a primitive lattice transla- 
tion. The existence of edge states is due to the combined 
symmetry. Since the edges of strip SSL2 are ferromag- 
netic and break the combined symmetry, the edge states 
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are gapped. However because the combined symmetry 
on the edges of strips CSL and SSL1 is still preserved, 
the edge states exist. But these edge states are no longer 
robust to disorder, which can be shown from conductivity 
calculations ((c) and (f) of Fig. ([3])). Here it is disorder 
that breaks the combined symmetry and the QSH phase 
in the system will no longer be protected. While for sys- 
tem with s z channel term, though ferromagnetic edges or 
disorder breaks the combined symmetry, the spin s z con- 
servation follows to assure the existence and robustness 
of the edge states. 

In conclusion, we introduce non-magnetic and mag- 
netic staggered potentials to a trivial insulator with spin- 
orbit coupling and find that they can induce a topologi- 
cal phase transition and drive the system into topological 
phase. For non-magnetic staggered potential, the result- 
ing QSH phase is protected by T symmetry and sup- 
ports edge state on any edge, which is robust to disorder. 
While for magnetic staggered potential, there is a com- 
bined symmetry of T and a primitive lattice translation 
in the system. If there is also an additional symmetry 
i.e., spin s z conservation, edge states will exist on any 
edge in band-inverted phase and are robust to disorder. 
However in the absence of such symmetry, though edge 
states will exist on specific edges in band-inverted phase, 
they are no longer robust to disorder. Our these results 
imply that though generally QSH effect is protected by 
T symmetry but if there are additional symmetries, QSH 
effect can also be found in T breaking systems. . 
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